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^ 1 INTRODUCTION 
(N 

^ Accreting stars with strong magnetic fields are generally observed 
I . to rotate more slowly than their less-magnetic or discless counter- 
y—^ ' parts. In protostars, T Tauri systems (which often have strong sur- 
IL" face magnetic fields of ~ 10^ — 10'^ G) with discs rotat e more 
. rapidly than systems without discs dCetman et al.l |2008|) . Most 
(but not all) mainsequence Ap stars (with surface fields of up to 
are observed to rotate much slower than normal A stars 
' ('Stepieii & Landstreetll2002h . and recent work has suggests this re- 
lationship extends d own to their pre-main-sequence progenitors, 
the Herbig Ae stars jAlecian et al. In high energy systems 

the result is similar: accreting neutron stars with weak (~ 10* G) 
fields rotate up to 10^ times faster than neutron stars with strong 
(~ 10^^ G) fields. These observations suggest that the interaction 
between the accretion disc and stellar magnetic field plays a critical 
role in regulating the spin-rate of the star 

Early theoretical studies of accretion predicted that a strong 
stellar field would truncate the accretion disc some distance from 
the surface of the star, with the truncation radius located roughly 
where the magnetic pressure {B^ /4-k) equals the ram pressure of 
the infalling gas (miir/Syrr^), so that infalling matter is channelled 
ont o the surface via magn etic field lines, causing the star to spin 
up. jPringle & Rees|[r972h . This assumes that the disc is truncated 
inside the corotation radius (rc = (GA/*/r2j) ^''''), where the 
star's spin frequency is equal to the disc's Keplerian frequency). 
If instead the magnetic field spins faster than the inner edge of 
the disc, a centrifugal barrier prevents accretion. Interaction be- 
tween the magnetic field and the disc will then spin down the star 



ABSTRACT 

We investigate the evolution of a thin viscous disc surrounding magnetic star, including the 
spindown of the star by the magnetic torques it exerts on the disc. The transition from an 
accreting to a non-accreting state, and the change of the magnetic torque across the corota- 
tion radius rc are included in a generic way, the widths of the transition taken in the range 
suggested by numerical simulations. In addition to the standard accreting state, two more are 
found. An accreting state can develop into a 'dead' disc state, with inner edge r-m well outside 
corotation. More often, a 'trapped' state develops, in which rin stays close to corotation even 
at very low accretion rates. The long-term evolution of these two states is different. In the dead 
state the star spins down incompletely, retaining much of its initial spin. In the trapped state 
the star asymptotically can spin down to arbitarily low rates, its angular momentum trans- 
feiTed to the disc. We identify these outcomes with respectively the rapidly rotating and the 
very slowly rotating classes of Ap stars and magnetic white dwarfs. 



MHD - stars: oscillations - stars: 



jlllarionov & Sunvaev 1975: Mi neshige et alJl991 : lLovelace et al.l 
ll999l : lRomanova et aU2004i : iUstvugova et alj2006l) . 

The presence of the centrifugal barrier is often equated in the 
literature with the idea that the accreting gas will be flung out, or 
'propellered' out of the system so as to maintain a steady state. 
This assumption turns out to be both arbitrary and unnecessary. For 
example, in order for the accreting material to be flung out of the 
system, the disc must be truncated a sufficient distance away from 
rc. Otherwise the rotational velocit y difference between the disc 
and the magnetosphere is too small dSpruit & Taamll9^3h . 

Steady disc solutions with a centrifugal barrier at t he in- 
ner edge were first described by ISunvaev & Shakural (ll977 V who 
called them 'dead discs', because even though the disc is actively 
transporting angular momentum outwards, no accretion onto the 
star takes place and the disc itself is very dim. 

As pointed out already in ISunvaev & Shakural ( Il977h and 
ISDruit&Taanil ( ll993l) . in a system with an externally imposed mass 
flux the likely effect of a centrifugal barrier is to cause the accretion 
onto the star to be cyclic. Accretion phases alternate with quies- 
cent periods during which mass piles up outside the barrier, without 
mass having to leave the system. In the quiescent phase, the angu- 
lar momentum extracted from the star by the disc-field interaction 
is carried outward through the disc by viscous stress. This alters the 
surface density profile of the disc from the usual accretin g solution. 

In our previous paper dP'Angelo & Spruid l2010 l hereafter 
DSIO) we studied this form of cyclic accretion with numerical so- 
lutions of the viscous diffusion equation for a thin disc subject to 
a magnetic torque. A s in the (somewhat more ad hoc) model of 
ISpruit & TaamI ( Il993h , limit cycles of the relaxation oscillator type 
were found. The cycle period of these oscillations depends on the 
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accretion rate, from fast oscillations at higher mass flux to arbitrar- 
ily long periods at low accretion rates. 

Instead of the two states: accreting and dead as suggested 
above, the results in DSIO are actually described better by includ- 
ing a third, intermediate state we call here the 'trapped' state: 

(i) rin<rc'- accreting state, star spins up, 

(ii) rc— A <rin<rc+A: trapped state, spinup or spindown, 

(iii) Tin-Vc^ A star spins down, no accretion (dead disc), 

where A <Crc is a narrow range around corotation, to be specified 
later. In state (ii), the inner edge of the disc remains close to corota- 
tion over a range of accretion rates onto the star, and the net torque 
on the star can be of either sign, depending on the precise location 
of the inner edge of the disc. 

For a given accretion rate, a disc that starts in state (i) will 
gradually move into state (ii) or (iii) as the star spins up and 
moves inward. In state (ii) accretion can proceed steadily or hap- 
pen in bursts, depending on the disc-field interaction at rin. For 
steady externally imposed accretion a disc in this state will eventu- 
ally move into spin equilibrium with the star, so that the net torque 
on the star is zero. In the dead state (iii) a steady state can exist if 
the torque exteited by the star is taken up at the outer edge of the 
disc by a companion star. If we neglect the transition to the pro- 
peller regime, then in theory the dead disc solution can exist for a 
disc truncated at any distance outside rc. Such a disc will remain 
static as the star spins down and rc moves outward. Our model 
is thus qualitatively different from the conventional 'propeller' pic- 
ture since at very low accretion rates a considerable amount of mass 
remains confined in the disc, and the star can be efficiently spun 
down. 

In the following we study the long-term evolution of the star- 
disc system by using the description of magnetospheric accretion 
in DSIO, allowing the star's spin rate to evolve. Of special interest 
will be the trapped state (ii), since in many cases the evolution of 
the system ends in it. The accretion cycles found in DS 10 also take 
place essentially within a trapped state. The inner edge of the disc 
is near corotation in the trapped state, as is the case also for a disc in 
spin equilibrium with the accreting star. Spin equilibrium is only a 
special case of a trapped state, however. In general a trapped state is 
not one of spin equilibrium, spinup is possible as well as spindown. 
We focus particular attention on cases in which the mean accretion 
rate through the disc is very low and the star cannot easily move 
into spin equilibrium with the disc. This is most applicable for iso- 
lated stars, or accreting binaries in quiescence. 

This scenario poses a number of questions which we address 
in the course of this paper. These include: under what conditions 
does the disc get into a trapped state, and when does it instead 
evolve into a dead state? It will turn out that this is determined 
by the details of the disc-field interaction and the ratio of the spin- 
down timescale of the star (Tsd) to the viscous timescale of the 
disc (Tvisc). The initial conditions of the disc also significantly in- 
fluence the outcome. In section [5!2l we ask how an initially trapped 
disc could become untrapped as a dead disc state. In particular, does 
this depend on the initial location of the inner edge of the disc, 
the initial accretion rate, the presence or absence of a companion 
or the size of the disc? Finally, in sec. |5] we discuss the physics 
that determines whether a disc will become trapped, and 15.3. ll we 
ask whether a trapped disc could plausibly regulate the slow spins 
observed in Ap stars, some of which have spin periods of up to 
decades. 

In a companion paper we investigate the observable conse- 
quences of a trapped disc, focusing in particular on how the burst 



instability studied in DSIO will change the spin evolution and ob- 
servable properties of the star. In that paper we also discuss our 
model's predictions in terms of observations of magnetospherically 
regulated accretion in both protostars and X-ray binaries. 

We use the code developed in DSIO, adding the star's moment 
of inertia as a parameter of the problem in order to follow the spin 
evolution of the star in response to the disc interaction. We can 
then simultaneously follow the viscous evolution of the disc and 
spin evolution of the star as the star's spin changes, and explore 
how these two interact with each other. We describe our model in 
more detail in the following section. 



2 MAGNETOSPHERIC INTERACTIONS WITH A THIN 
DISC 

2.1 Magnetic torque 

The interaction between a strong stellar magnetic field and sur- 
rounding accretion disc truncates the disc close to the star, and 
forces incoming matter to accrete along closed field lines onto 
the surface of the star in a region called the magnetosphere. At 
the outer edge of the magnetosphere (termed here the magneto - 
spheric radius), the field lines become strongly embedded in the 
disc over some small radial extent that we term the interaction re- 
gion, Ar. The differential rotation between the star and the Keple- 
rian disc will cause the field lines to be twisted, which will generate 
a toro idal component to an initially poloidal field (e.g. lGhosh et al.l 
Il977h . This will allow the transfer of angular momentum between 
the disc and star, with the torque per unit area exerted by the field 
on the disc given by r = rSz^z, where: 

Sz4' = (1) 

is the magnetic stress ge nerated by t he twisted field lines. Both the- 
oretical arguments (e.g. Iai V 1985; 'Lovelace et al. 111995*) and nu- 
merical simulations (such as Miller & Stone 1997; Goodso n et al.l 
ll997l ; lHavashietal.ll996l) suggest that the strong coupling between 
magnetic field lines and the disc will cause the field lines to inflate 
and open. The inflation and opening of field lines limits the growth 
of the B^ component for the field to B^ — rjB^ , with rj of order 
unity, and reduces the radial extent of the interaction region, since 
beyond a given radius the field lines are always open and the disc- 
field connection will be severed. We take the interaction region to 
be narrow, Ar/r < 1 (as found in numerical simulations, see sec- 
tion 2 of DS 10 for a more detailed discussion). Assuming the star's 
dipole field strength Bd {r) as an estimate of Bz, and taking into ac- 
count that S acts on both sides of the disc, l[T) yields the magnetic 
torque To exterted on the disc: 

To = 47rrArrSz</, = ryr^Arfii (2) 

This torque exists only if the inner edge rin of the disc is out- 
side the corotation radius rc. For rin < rc, we have instead a disc 
accreting on an object rotating slower than the Kepler rate at rin. 
By the standard theory of thin viscous discs, the torque exterted on 
the disc by the accreting object then vanishes, independent of the 
nature of the object. The torque Tb (ri) thus changes over a narrow 
range around rc. To model this transition we introduce a 'connect- 
ing function' j/e : 

Ts(rin) =j/E(nn)ro(ri,), (3) 
with the properties j/s ^ (rc — rin 2> Ar), ^ 1 (rin — rc 2> 
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Ar). As in DS 10, we take for this function 



j/E = — 1 + tanh 



Ar 



(4) 



Tlie widtli of the transition is thus described by Ar. We take the 
same value for it as used in eq. (|2j- 



2.2 Model for disc-magnetosphere interaction 

In DSIO we derived a description of the interaction between a 
disc and magnetic-field for a disc truncated either inside or out- 
side rc, and introduced two numerical parameters to connect the 
two regimes. To keep the problem axisymmetric, we assumed a 
dipolar magnetic field, with the dipole axis aligned with the stel- 
lar and disc rotation axis. Since the region of interaction between 
the disc and the field is small, we use our description of the inter- 
action as a boundary cond ition for a standard thin accretion disc 
dShakura & Sunvaevill973h . 

To evolve a thin disc in time, we must choose a description 
for the effective viscosity [v) that allows transport of angular mo- 
mentum. We adopt an a prescription for the viscosity and assume 
a constant scale height (/i) for the disc, so that: 



v^a{GM,f^''{h/rfr^'^. 



(5) 



At the inner edge of the disc the behaviour is regulated by the disc- 
field interaction. However, since the interaction region is small, we 
incorporate the interaction as a boundary condition on the inner 
disc, and assume that the majority of the disc is shielded from the 
magnetic field and then evolves as a standard viscous disc, albeit 
with a very different inner boundary condition from the standard 
one. Below we summarize our analysis of the disc-field interaction 
and how these translate into boundary conditions on the disc. (For 
the detailed derivation of our boundary conditions, see sections 2.3, 
2.4, and 3.2 of DSIO). 

2.2. 1 Surface density at rin 

In a dead disc, the disc-field interaction prevents matter from ac- 
creting or being expelled from the system, instead retaining matter 
that interacts with the magnetic field. This implies that the angular 
momentum injected via magnetic torques in the interaction region 
Ar must be transported outwards by viscous torques in the disc. A 
dead disc will therefore have a maximum in surface density at rin, 
and S(rin) will depend on the amount of angular momentum being 
added by the disc-field interaction. 

By equating the amount of angular momentum added by the 
field to the amount carried outwards by viscous processes, we 
can calculate the surface density at the inner boundary of the disc 
needed for the injected angular momentum to be carried away. This 
yields (see DS 10) a value for the surface density E at the inner edge 
of a dead disc, proportional to the magnetic torque (|3): 



37ri/E(rin) 



Tb 



(6) 



where JIk is the Keplerian rotation frequency. If the stellar field is 
a dipole and we use (|5]l to describe the viscosity, then for rin > rc, 



E(ri.O 



. E(rin) thus decreases rapidly with increasing rin. 



2.2.2 Accretion rate across rc 

If the inner edge is well inside the corotation radius rc, we use a 
standard result to estimate the location of rin as a function of the 



accretion rate m^. It is obtained from the azimuthal equation of 
motion for gas at the point at which it is forced to corotate with the 
star (c.f. lSpruit & Taamll 19931) . This gives: 



rin7r(5'z0)/r2, = m. 



(7) 



[Note that we take the sign of rh positive for inward mass flow.] It is 
not necessary that stationarity holds: (|7} can also be applied when 
the inner edge of the disc moves. However, since it describes the 
accretion through the magnetosphere-disc boundary rin, it has to 
be applied in a frame comoving with rin. If m is the mass flow rate 
in a fixed frame, it is related to the accretion rate in this comoving 
frame (co) by 



= m + 27rrinE(rin) rii, 



(8) 



where fin is the rate of change of the inner disc edge. 

To connect the accreting case with the dead disc case, for 
which m = 0, we need one more prescription, this time for the 
accretion rate as a function of the inner edge radius. We introduce 
a connecting function ym for this (DSIO): 



with the properties y, 
Ar2), with 



mco(rin) = J/m(rin)rha(rin), 

1 (rc-rin > Ar2), ym 



1 — tanh 



Ar2 



(9) 

0(rin-rc > 
(10) 



where Ar2 describes the width of the transition (different in general 
from Ar). 

With the star's assumed field of dipole moment /i, = ^t/r^ 
and Keplerian orbits in the disc, l[8j becomes, with the viscous thin- 
disc expression for rh: 



„ \ 2 d 1/2 

ar 



)L =?/" 



27rrinE(rin)fin, (11) 



where r; is a numerical factor of order unity and /i the dipole mo- 
ment of the star (see DS 10 for details). 

Along with our description for the viscosity, (|6) and (II U de- 
fine a boundary condition at rin and an equation for rin(t), for a disc 
over a continuous range of accretion rates, from strongly accreting 
systems (rin <S rc) to dead-disc systems (m ~ 0). 

2.2.3 Evolution of corotation radius 

In order to study the response of a disc to changes in spin of the star, 
we must incorporate the angular momentum exchange between the 
star and disc: 

dn« _ dJ 

At ^ At' 

which introduces the moment of inertia of the star, /* = k^A{t.R^ 
as an additional parameter of the problem. 

The disc-star angular momentum exchange AJ/At has two 
components: matter accreting onto the star adds angular momen- 
tum at a rate rh,corfni^K(rin), while the disc-field coupling outside 
co-rotation extracts angular momentum spinning the star down. 
The rate of angular momentum exchange between the disc to the 
star will thus be (with|2l|3]|4j: 

dJ 
At 



L- 



(12) 



mcorinriK(rin) - Tb 



^3 



3/2 



Ar 



(13) 
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The corotation radius (a function of fi*) evolves as: 



At 



2dJ 1 f GM, 
' 3 * V ^5 



-1/2 



(14) 



Eq. U3t shows that there is a value of rh for which there 
is no net angular momentum exchange with the star. This is the 
'spin equilibrium' state discussed in previous work. This zero-point 
will depend on our adopted connecting functions, as well as the 
size of the transition widths, Ar and Ar2. If m = 0, there is no 
spin-equilibrium solution: the star will spin down by the magnetic 
torque. 

2.2.4 Steady-state solutions 

In the presence of a magnetic torque at the disc inner edge, the 
steady solutions (d/dt = 0) of the thin viscous disc diffusion equa- 
tion with the above boundary conditions have the form (cf. DS 10): 



Tb 



nn\l/2 



+ rri 



1 - 



1/2 



(15) 



where m is the accretion rate onto the star, given by (|9). If the 
inner edge is inside corotation {Tb = 0), S has the standard form 
for steady accretion on an object rotating below the Keplerian rate 
(second term on the RHS). 

For Tin well outside corotation (vin-rc^ Ar), m 4- and we 
have a dead disc. The surface density is then determined by the first 
term on the RHS. The steady outward flux of angular momentum in 
this case has to be taken up by a sink at some larger distance, oth- 
erwise the disc can not be stationary as assumed. This sink can be 
the orbital angular momentum of a companion star, or the disc can 
be approximated as infinite. The latter is a good approximation for 
changes in the inner regions of the disc, if time scales short com- 
pared with the viscous evolution of the outer disc are considered. 

2.3 Numerical method 

We use the one-dimensional numerical code described in DSIO to 
evolve the standard diffusive thin-disc equation with our viscos- 
ity prescription (|5) and our description of the disc-field interaction 
(which gives the inner boundaiy conditions the boundary condi- 
tionsl6land lI U . At the outer boundary a mass flux and a flux of an- 
gular momentum are specified in various combinations (described 
in sec. 12.3. iT ). 

The calculations are done in dimensionless coordinates and 
variables. In DSIO we scaled all physical lengthscales to rc, and 
physical time scales to Tvisc(?'c). Since in this paper we want to 
follow the evolution of r^, we instead use the stellar radius r* and 
Tvisc{r*) = Tt to scale our physical length and timescales. The 
grid is logarithmically spaced (to ensure sufficient resolution in the 
inner disc to capture the disc instability). It is an adaptive mesh, 
such that the inner boundary moves with rin. 

Since the grid used is time dependent, the outer boundary con- 
dition is also applied at a time-varying location. As discussed in 
DSIO, the artefacts this causes are small, compared to specification 
at a fixed location (at least for the large discs studied in most cases). 

The size of the discs studied range from 10 to 10® times the 
inner edge radius, and the number of grid points needed for suffi- 
cient resolution varies accordingly, from 90 for the smallest to 560 
for the largest discs. 

Time-stepping is done with an implicit method, so that the 
short timescales encountered during episodes of cyclic accretion 



can be followed, as can the much slower viscous evolution of the 
disc as a whole and the spin down the star. It is adapted to the stiff 
nature of the equation to be solved (see DS 10 for details). 



2.3.1 Outer Boundary Condition 

The lifetime and evolution of a star surrounded by a dead disc is an 
inherently time-dependent problem, so the initial conditions in the 
disc can be critical for its evolution. Since the spin-down timescale 
for the star can be much larger than viscous timescales throughout 
the disc, the conditions in the outer disc will also strongly influence 
the evolution of the system. 

We thus consider the effect of varying the outer boundary con- 
ditions for the disc. The first condition we study is the simplest: a 
fixed mass flux m = rho (> 0, corresponding to accretion). As dis- 
cussed in the introduction, a key aspect of disc-magnetosphere in- 
teraction is that accretion is possible even as the star is spun down. 
At fixed m > 0, the angular momentum flux can be either inward 
or outward. 

If the mass flux specified vanishes at rout, the boundary con- 
dition is 



dr 



(r'/'i/E)| 



0. 



(16) 



On long evolution timescales, the finite extent of the disc it- 
self could be relevant in a star without a companion where the disc 
can spread outwards. To model this, as our final boundary condi- 
tion we take S(rout) ~ 0, so that the angular momentum added at 
rin is carried away by the outer parts of the disc, causing the disc 
to spread outwards. In section l?!2l we discuss the consequences of 
these assumptions in limiting the lifetime of a trapped disc. 



2.3.2 The Evolution of rc 

The final modification to our code used in DSIO is to allow to 
evolve as the spin rate of the star changes JI4t . The characteristic 
evolution timescale for rc, the spin-down timescale for the star, is 
much longer than the nominal viscous timescale in the disc (see the 
next section). The code updates by an explicit time step, rather 
than implicitly, as we do the other variables. Rather than discretiz- 
ing ( I14t and add it to our system of linearized equations that are 
solved numerically at each timestep, we instead approximate the 
evolution in rc to first order in time, that is: 



rc(fo + At) = rc(to) + 



drc 
'At 



At 



(17) 



to 



This scheme is simpler than adding additional equations to the 
code, and is sufficient to describe the co-evolution of the disc and 
stellar spin-rate. However, as we discuss in sec. 14.31 it is not ac- 
curate enough when rin is close to rc and the spindown timescale 
comparable to the viscous timescale in the inner regions of the disc. 



3 CHARACTERISTIC TIMESCALES OF DISC-STAR 
EVOLUTION 

Three kinds of timescale play a role in the evolution of a star cou- 
pled magnetically to an accretion disc. These are the spin period of 
the star, the timescale for changes in spin period of the star, and vis- 
cous evolution timescales of the disc. The viscous evolution does 
not have a single characteristic timescale; it can vary over many 
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orders of magnitude depending on which regions in the disc partic- 
ipate in the evolution. 

The spin period of the star (P,) determines the location of the 
corotation radius which in turn sets the accretion rate at which the 
magnetic torque changes sign. P, also determines the timescale 
for magnetospheric variability (from processes like reconnection 
of field lines), which can lead to variability in m, Ar, Ar2 and 
the component of the magnetic field (which sets the magni- 
tude of the torque). P, is much shorter than the other characteristic 
timescales studied in this paper, and the complex variability pro- 
cesses are best studied with detailed MHD simulations, so in this 
work we assume time-averaged values for m, Ar, Ar2 and 
and neglect timescales much shorter than the characteristic viscous 
timescale at P, . 

A convenient unit of time for measuring changes in a viscous 
disc at a distance r from the center is ta = r"^ /v{r), sometimes 
called the accretion- or viscous timescale at distance r. If a is the 
viscosity parameter and H the disc thickness, it is longer than the 
orbital timescale Q,^ by a factor {r / H)'^ , a large number for 
most observed discs. Natural choices for r in this expression would 
be the disc's inner edge rin (tin) or the corotation radius rdtc)- 
Both of these are functions of time. The actual timescales of varia- 
tion in our discs can be much shorter than tc and U-n, however, since 
the extent of the disc that participates in the variation can be much 
smaller than rin. In the cyclic accretion mode described in DSIO, 
for instance, cycle times as short as O.Olic are found. The timescale 
for viscous adjustment in the outer disc regions, on the other hand, 
can be very large compared to tc- 

The longest timescale is the rate at which the star's spin will 
change, which is determined both by the rate of angular momentum 
exchange with the disc and the star's moment of inertia. The spin- 
down torque of a dead disc (with m = and rin = Vc) is, from eq. 
GD: 

= (18) 

At rf„ 

where 5 = Ar/rin. The characteristic spindown time is: 

TsB = p,/p^ ;7 . (19) 

If the inner edge stays near corotation, Qk = 5^*, replacing ri„ by 
Tc, this yields: 



quantity Tvisc/Tso thus defines the ratio of the time that gas at that 
radius takes to travel inwards onto the star and the time it would 
take to spin-down so that rc = rin, independent of radius. In all 
cases, the spin-down timescale is much longer than the viscous 
timescales in the disc, so that at least part of the disc is able to adjust 
to the new spin rate of the star. However, the exact ratio between the 
two timescales varies over ten orders of magnitude, from 0.06 for 
a strongly magnetic Herbig Ae star to 10"^^ for an accreting mil- 
lisecond pulsar This ratio implies that the extent of spin-down will 
be influenced by the viscous evolution of the disc itself in response 
to the disc-field interaction, and that this evolution will vary sub- 
stantially between different systems, breaking the scale invariance 
usually assumed in disc-magnetospheric interactions. In section l431 
we demonstrate how the ratio of these two timescales is critical in 
determining the ratio of rin to Vc in a trapped disc. 

3.1 Representative Model 

In sections |4] and 15.21 we study how a trapped disc can form and 
evolve, as well as how it can become untrapped. In order to sim- 
plify comparison between different simulations, we adopt a set of 
parameters for a representative model, which we then vary between 
solutions as necessary. 

For the dimensionless parameters we adopt Ar/rin = 0.1, 
Ar2/rin = 0.04, and P^/B^ = V = 0.1. The values of Ar/n^ 
and Ar2 / rin are small enough to provide an abrupt transition be- 
tween an accreting and non-accreting disc, but do not show the 
cyclic instability discussed in DSIO. Neglecting the star's spin 
change, the problem has a scale invariance (DSIO), whereby the 
parameters ^, Af«, fl* and m can be re-written as the ratio rh/rhc, 
and The is the accretion rate in (|7) that puts the magnetospheric ra- 
dius at Tc, a natural unit of m for magnetospheric accretion. The 
results in DS 10 were presented in this unit. 

As discussed in the previous section, the variation of r-cwith 
time during spindown of the star makes this unit impractical. In- 
stead, we present the representative model in units suitable for a 
protostellar system with Tvisc/Psd = 2.6 x 10"'' (which is as large 
a ratio as the present version of the code allows), and explore the 
effect of varying this ratio. As unit of length we use r*, the star's 
radius, and for timescale the nominal viscous timescale of the 
disc at the star's radius. 



^- = 1^^" ("K(nn) = f^,). (20) 

The spindown timescale varies considerably between different 
sources. Adopting = 1 and Ar/rin = 0.3, this spindown 
timescale is short enough to account for spin-regulation in slowly 
rotating magnetic stars. In Table |3] we summarize the predicted 
spin-down timescales for a slowly rotating X-ray pulsar, a millisec- 
ond pulsar, a slowly rotating Ap star, and a typical T Tauri star. 
For all these examples but the millisecond pulsar the spin-down 
timescale is much shorter than the lifetime of the star. Provided the 
conditions are such that the inner edge of the disc can stay near 
corotation (i.e. what we have called the 'trapped disc' state), it will 
be able to spin down a star to very slow rotation periods. In the next 
sections we will explore how this could work in detail by evolving 
a viscous disc in time numerically. 

The last column of Table [3] lists the ratio of the viscous 
timescale r^/u and il9i . Note that for our description of viscos- 
ity, the viscous and spin-down timescales both scale as r"^''^. The 



4 TRAPPED DISCS 

4.1 Trapped disc evolving from an accreting disc 

For our description of the disc-field interaction (which ignores 
outflows), once the accretion rate falls to zero, the inner edge of 
the disc could be located anywhere outside rc, depending on the 
amount of mass in the disc. What then determines the location of 
the inner radius of a dead disc? To answer this question, we simu- 
late an initially steadily accreting disc in which the accretion rate 
at rout suddenly decreases to zero. As the reservoir of gas in the 
disc runs out, the accretion rate onto the star declines, and the inner 
radius of the disc moves outwards. 

In the simulation we use our representative disc parameters 
described above, and set the initial inner radius of the disc to be 
just inside rc, rin(t = 0) = 0.88rc,o, and the outer radius rout ~ 
100 rin. We can calculate the corresponding accretion rate from (|7), 
and use the static solution for E given by l |15l l as our initial surface 
density profile. At t = 0, we set 9r(i^Er^^^)|rout = 0, so that no 
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Table 1. Spindown and viscous timescales for different type of magnetic stars 



Source 


Mass 


Radius 


-B* 


P* 






T'visc(»'c)/TsD 




(Mq) 


(R©) 


(G) 






(years) 




slow Pulsar 


1.4 


1.4 X IQ-"' 


1012 


5.0s 


2.9 X 10-" 


4400 


3 X lO-'^ 


ms Pulsar 


1.4 


1.4 X 10-5 


10* 


0.1s 


2.9 X 10-" 


2.6 X 10" 


2 X 10-1^ 


Magnetic Ae star*^ 


3.0 


5.5 


10* 


10 yrs 


4.0 


3 X 10^ 


0.06 


T Tauri Star_^ 


0.6 


3.0 


1500 


7 days 


0.54 


2.3 X 10"' 


0.001 



" B, and I, fromlStepie^ j200(ll) 
* ISiposetaUiiool 



mass is added to the disc or allowed to escape. This sets a constant 
angular momentum flux at rout. 

The results are shown in Fig. [T] The bottom panel of Fig. [T| 
shows the change in accretion rate onto the star, scaled to riic. The 
top panel of Fig. [T] shows the evolution of the inner radius (solid 
black curve) and (dashed red curve) in response to the changing 
accretion rate. After initial steady accretion over about 30 t* (less 
than 1/lOth the viscous timescale at rout), rh through the inner edge 
of the disc begins to decrease as the reservoir of mass in the disc 
is accreted onto the star, and rin begins to move outward. From 
30 — 1500 T* , m decreases exponentially with a decay timescale of 
about 240 t, and the disc moves outwards. However, rin increases 
by only about 20% as the accretion rate decreases by three orders of 
magnitude. (The structure of rin around rc is an artifact of the tanh 
connecting functions we adopted to describe rin and rh across the 
transition region). 

After ~ 1500 t, the star begins to spin down (so that rc moves 
outwards), and the inner radius of the disc begins to track rc, so that 
the ratio rin/rc remains nearly constant thereafter. The behaviour 
of the accretion rate at this point also changes. Although it contin- 
ues to decrease exponentially, the decay timescale lengthens con- 
siderably and the accretion rate (~ 10~* of the initial rh) is reg- 
ulated by the spin-down rate of the star. Instead of continuing to 
move away from rc into the 'dead disc' regime (in which rh = 0), 
the inner radius instead remains trapped at nearly a constant frac- 
tion of rc while the star continues to spin down. We thus call this 
disc solution a 'trapped disc", since rather than continue to move 
outwards, rin becomes trapped at a nearly constant fraction of rc. 

At the outer edge of the disc m = and there is an outward 
angular momentum flux. The accretion onto the star comes from 
the disc being slowly eroded (although at a very low rate) as rin 
moves outward. The evolution of rin and the inner parts of the disc 
is determined by the spin down rate of the star itself, which is itself 
influenced by how close rin can stay near rc. In Sec. l4.3l we demon- 
strate how rin/rc is mostly determined by the parameters Ar and 
Ar2, and the ratio TvIsc/Tsd- However, the main conclusion of 
this section is clear: if a trapped disc forms and can efficiently carry 
away the angular momentum of the star, over spin-down timescales 
the disc will accrete at such a rate so as the inner edge of the disc 
can move outwards together with rc and the star could in theory 
spin down completely. 



4.2 Trapped disc evolving from a dead disc 

Consider next a case where the initial condition is a dead disc, 
(rin > rc) given by the steady profile ( dlSI l with rh = and 
= !)• As the star spins down, rc moves out until it catches up 




10^ 10^ 
Time [T^i, 



Figure 1. Transition from an accreting to a 'trapped disc' state. The initial 
E profile of the disc sets rin = l-5r«, but the accretion rate through rout 
is set to zero at t = 0. As a result of mass loss by accretion through rc, 
the disc quickly evolves away from a steady accretion state and after about 
1500 tt settles into a slowly evolving state in which rin tracks rc. Top: 
The black solid curve shows the evolution of the inner radius in time, while 
the red dashed curve shows the evolution of rc. Bottom: The accretion rate 
onto the star, which decreases sharply as rin moves outwards across the 
corotation radius. 



with the inner edge rin. From then on, the same evolution is as in 
the previous case: rin and rc move outward together indefinitely. A 
small amount of mass accretes while the star's angular momentum 
is transferred to the disc. 

The asymptotic evolution of this dead disc can be compared 
with the case when a fixed mass flux is imposed at the outer edge. 
The asymptotic state is then a steady state with spin equilibrium: 
the spin-up torque due to the accreted mass is balanced by the mag- 
netic torque at rc transfering angular momentum outward. 

We illustrate the distinction between these two cases in figure 
[2] This shows the evolution of an initially dead disc (black, thick 
curves) and accreting disc (red, thin curves). Both discs have the 
same representative parameters (Sec. 13. it and rout = 100 rin, but 
with different initial inner radii (a few times the stellar radius), ac- 
cretion rates and appropriate initial surface density profiles given 
by dlSb . For the dead disc, we take rin = l-3rc,o (where rc.o 
is the initial corotation radius) which corresponds to rh ~ for 
our chosen value of Ar2/rin. For the accreting disc, we choose 
Tin = l.lrc.o, which corresponds to a low but non-zero accretion 
rate (rh = 8 x 10"^rhc). 

Since rh for the accreting disc is initially low compared to 



Long-term evolution of discs around magnetic stars 7 



DC 



DC 
DC 



2 
1.4 
1.3 
1.2 
1.1 
1.0 




10" 



10' 
Time [J,, J 



10" 



10' 



Figure 2. Comparison of the evolution of rin and rc between a dead disc 
(black solid and dashed curves) and accreting disc (red dot-dashed curves). 
Top: The evolution of rc and rjn in units of the stellar radius. Black curves: 
Tin (solid) and rc (dashed) in the dead disc; red curves: the accreting disc 
(rin, dot-dashed curve; rc, triple-dot dashed curve). Bottom: The ratio 
I'ln/fc for the dead disc (solid) and the accreting disc (dot-dashed). The 
dead disc keeps evolving indefinitely, the accreting case reaches a steady 
state in spin equilibrium with the star around t ~ 10^. 



rhc, at early times rc evolves at the same rate in both simulations 
(accreting: red, triple-dashed curve; non-accreting: black, dashed 
curve), and the star spins down. Eventually, however, the amount 
of angular momentum added by the accreted gas becomes compa- 
rable to the amount removed, and spin equilibrium is reached at 
~ 10^ tt ■ In contrast, for the m = case, the disc at first remains 
unaffected, while rc moves outwards. (This is because the magnetic 
torque depends only on distance, not on the rotation rate of the star). 
After Tc moves close enough to rin that accretion can begin (around 
4 X 10'* t*), the two start to move outwards at approximately the 
same rate. The (low) accretion rate onto the star is determined by 
the (slow) rate at which rc moves outwards, and the star continues 
to spin down indefinitely. The bottom panel of fig.[2]shows the ratio 
of Vin/vc for the accreting (dashed) and non-accreting disc (solid). 
After the non-accreting disc passes out of the dead disc phase (at 
~ 4 X 10^), in both systems the ratio changes by less than 10%, 
and Tin always remains close to rc. The non-accreting disc, how- 
ever, differs in that it never reaches spin equilibrium. 

The main difference between the evolution of a dead disc and 
an accreting disc is the behaviour of the inner edge radius. As seen 
in sec. 14. 21 in the initially dead disc the accretion rate onto the star is 
determined by the disc's behaviour when it reaches rin — rc, while 
in an accreting disc the accretion rate is governed by the value set 
at the outer boundary. 

Both the accretion rate on the star and the outward angular 
momentum flux in our trapped discs depend sensitively on the dis- 
tance between rin and rc compared with the transition widths Ar 
and Ar2. In the results of figure|2] rin- rc is of the order 0.5 — 2 Ar. 
In the next section we develop an analytic estimate of this number 
and compare it with the numerical results. 



4.3 Analytic estimates for a trapped disc 

As we showed above, an initially dead disc will eventually start ac- 
creting at a low rate, in such a way that rin moves outwards together 
with rc at a nearly constant ratio. The accretion rate onto the star is 



determined by how close rc can move to rin before the disc moves 
outwards in response. The actual distance on which rin settles in 
cases like those show in the previous section depends on the details 
of the disc-field interaction (namely the parameters Ar and Ar2) 
and the ratio of timescales, Tvisc/TsD. 

The spin-down timescale derived above assumes that rin 
moves steadily outwards at the same rate as rc. This timescale is 
an upper limit, since as rc approaches rin there is reduced trans- 
port of angular momentum through rin and accretion onto the star 
begins. In addition, in order for spin-down to remain efficient, the 
angular momentum added by the disc-field interaction can be trans- 
ported through the disc and carried away at rout, otherwise rin will 
move quickly away from rc and spin-down will effectively cease. 

When rin is far enough from rc that m ~ 0, rin stays fixed 
as the star is spun down and rc moves outwards ( I15t . However, 
once rin moves closer to rin (within Ar or Ar2), this static state 
is no longer possible: either matter at rin starts accreting onto the 
star (ym 7^ 0), or the surface density at rin declines (j/e 7^ 1), 
which causes rin to move closer to rc until accretion through rin 
can begin. Since the viscous timescale in the inner part of the disc is 
much shorter than the spindown timescale, after accretion through 
rin begins, a pseudo-steady-state develops, and rin moves slowly 
outwards with rc. If the disc can maintain a steady-state for the 
given m, then rin will track rc, and the disc will remain a nearly 
dead disc as the star spins down to a small fraction of its initial spin 
period. 

We can study this quantitatively by considering the equations 
for rin and rc. The evolution of the inner edge of the disc dl llfrom 
Sec. 12. 2. 2t . defining it = Er for convenience is: 



27rM(rin)rin = Vm 



4r2*rf, 



du I 
dr ' 



(21) 



As long as the two terms on the right of M\\ balance, rin = 
even after accretion through rin begins. This will continue until the 
surface density profile near rin no longer satisfies dlSt . Since the 
change in rc is the only source of variability in the problem, rin 
will approximately track rc. 

Eq. l l21t cannot be solved as is, since du/dr depends on solv- 
ing the full time dependent diffusion problem. As an estimate we 
assume that u changes as a result of the changing boundary con- 
dition (which will increase the surface density gradient), divided 
by the rate at which the rest of the disc can respond to that change 
(which will smooth it out). This can be approximated by: 



du 

dr ' 



du{rin) _i 

V ■ 

dt 



(22) 



where v^isc is the viscous speed at rin, of order: v^isc ~ T^/r. 

The time derivative for ti(rin) follows from the boundary con- 
dition for u: 



du I 



du . du . 

■7: ^in "1" Ti '"c 

arin drc 



(23) 



The equation for rin then becomes: 



Tin = 



i^oymf^^ fAr 
inndn^yT, 



-1/2 

r.., — r. 



10 



3 dn^ysr in 



- 1 



(24) 

where / = (rc/rin), and we have used the definition of ys from 
Sec. 12.2. fi so that dr^^vr, = —dr^ys- 

The evolution of rc depends on the rate of angular momentum 
exchange with the star. Matter falling onto the star spins it up, while 
the interaction with the disc outside rc transfers angular momentum 
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outward and spins the disc down. The equation for this evolution is 
given by J14t . which can be re-written: 



2 



3 (GAf.)i/2/, 



7/2 



Ar ym rl/2 

— - -r-f 

rin 4 



'1/2 



(25) 



We can study the evolution of rc and fin in two limiting cases. 
In the limit where rin — rc 3> Ar, Ar2: 



Ar^7/2 



so that: 

rin = 



3 {GM,y/^h rin " 
0, 



(26) 



^w'Zin^Ar 5/2 
' c,0 



3 (GM*)i/2j, 



-2/5 



(27) 



rin.O- 



This is the limiting 'dead disc' case, where the amount of an- 
gular momentum being injected at rin can be extracted at rout and 
the disc remains steady while the star is spun down. It predicts a 
slightly smaller spin-down torque than was estimated in Sec.[3]be- 
cause the torque scales with rin. 

The inner radius of the disc will remain approximately con- 

., . , Vm drVS I Ar 

stant until either - — or — — ' 



become non-negligible 

J/E J/E V ''in y 

( I24t . Based on our assumption that the disc will remain in a quasi- 
steady-state while Tin moves, rin will evolve only in response to 
changes in rc, which means that / is a constant. This simplifies 
( 124b and ( 125b considerably: 



rin = 

where. 

Bo = 
Bi = 



(28) 



3 {GM^y/^L rin 



(29) 



-y-n 



6 (GAf.)i/2/, 

uoym ( Ar 



8rin9ri„(yE) \ rm 



3rinft i„(yE) ' 

The solution is then: 



(I {A,f^A,fl-)t 



3/2 



2/3 



(30) 



rin = 



3 (B of^^ -A^fl^ + Arf 
Bi-1 



t + r 



3/2 



2/3 



We can use ( |28b to calculate /, that is, how close rc can move 
towards rin before the disc will start moving outwards in response. 
Setting fc = /fin, we can re-express the constants in ( 130b as: 



Ar 
2 — /J/E 

' in 



lOj/E 



3rin9ri„yE 
rin9ri„ys 

ym 



+ ,/-! 
3 TsD 



(31) 



STvi 



and solve / numerically to give the approximate evolution for rc 
and rin in time. 

In Fig. [3] we compare our estimates for the evolution of rc 
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2.5 - 



2.0 - 



1.5 - 
1.35 



1.30 - 




10^ lO'^ 
Time [T, J 

Figure 3. Comparison between numerical solution and analytic estimate 
for an idealized dead disc around a star spinning close to break-up rc= 1.7, 
with Ar/rin = 0.1 and Ar2/rin = 0.05. Top: Numerical solution of the 
evolution of rc (black, dashed curve) and Tin (black, solid curve) in time. 
Oveiplotted is our analytic estimate for rin (red dash-dotted curve) and rc 
(red dash-triple dotted curve) for the same physical conditions. Bottom: The 
ratio rc/rin for our numerical solution (solid curve) and analytic estimate 
(dash-dotted curve). The disagreement in the two solution arises from our 
simplified treatment of the viscous disc. 



and rin to the solution from our numerical simulation. We con- 
sider a rapidly spinning star (rc = 1.7), with Ar/rin ~ 0.1 and 
Ar2 / rin = 0.05. In the top panel we compare the two solutions for 
rin (estimate: dot-dashed red curve vs. numerical solution: black 
solid curve) and rc (estimate: triple-dotted-dashed red curve vs. nu- 
merical solution: black dashed curve). In the bottom panel we plot 
the ratio rin/rc for the numerical (solid curve) and analytic (dashed 
curve) result. At early times, rc < rin ( 127b . and the solutions match 
exactly. However, at late times the solutions disagree somewhat. 
Most obviously, the value for / calculated by ( 130b is smaller than 
the numerical solution, that is that fin = for longer than we pre- 
dict, and there is some evolution in / over long timescales. This 
mismatch comes from our simplified treatment for 9r?i|^ , which 
over long timescales will depend on the surface density gradients 
in the entire disc. 

By comparing the sizes of each term on the left-hand size of 
( 131b . some insightful approximations can be made. Since the solu- 
tion is nearly a dead disc, the disc's accretion rate will be very low, 
so that 2/m <C 1, J/E ~ and dn^ys <C 1. As well, Hn is close 

to rc so that / ~ 0(1). Since in general Tvisc/TsD ^ 1, Olb can 



be approximated: 



160 Ar^yl 

9 ^in ym 



(32) 



and used to estimate /. The left-hand side of J32b is dominated by 
y,7/, which quickly grows as rin moves away from rc, and must 
balance the right-hand side. The ability to sustain a dead disc will 
thus depend on both the interaction between the disc and the field 
(through the parameter Ar2), and the ratio between the spin-down 
timescale and the accretion timescale. 

For a given TsD/Tvisc, if Ar2/rin <C 1, then ym is nearly a 
step function, and / will stay close to 1 even if TsD/2visc is very 
large. On the other hand, if Ar2 /rin 1, and accretion continues 
even when rin moves a fair distance from rc, then it is possible 
that the solution to ( 131b will predict a larger value for rin than can 
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support a dead disc. In this case no dead disc solution exists: even if 
tiie disc initially begins as a dead disc, once Vc moves close enough 
to Tin that accretion begins, the disc will begin to move outwards 
until matter at rin can be expelled from the system. The outflow of 
material will then proceed at a moderate rate as the surface density 
profile of the disc evolves away from the dead disc solution [ dlSt 
with m = 0] to the standard disc solution, with outflow rather than 
accretion onto the star 



5 TRAPPED AND UNTRAPPED 

5.1 Accreting discs evolving to trapped or dead disc states 

The results found so far and in DSIO indicate the strong tendency 
for the inner edge of the disc to track the corotation radius, what 
we call here a 'trapped disc'. This does not happen in all cases, 
however We would like to find out under what conditions a disc 
gets stuck in this way, and when instead the inner edge proceeds to 
move well outside corotation into the 'dead disc' state. 

Armed with qualitative understanding from the previous sec- 
tion, we can address this question with a few numerical exper- 
iments. We take the case of a neutron star with field strength 
Bs = 10^^ initial spin period P* = 5s, and initial inner edge 
radius rin = 0.95rc,o. The disc is thus initially in an accreting 
state. The (initial) outer boundary is located at rout ~ lOOrin, and 
we set m = there so that no matter can escape the system. The 
other parameters are the same as for our representative model. 

We first investigate the effect of varying the viscosity on the 
way the transition from an accreting to a dead disc takes place. 
This is shown in fig. ID The spindown timescale l |20| ( for the initial 
rc is TsD ~ 10^ years. From top to bottom, we plot this transition 
for decreasing values of viscosity (and hence increasing viscous 
timescales). To compare these two quantities we define Tvisc as the 
viscous timescale at the initial Tc. In the different plots, the ratio 
Tvisc/TsD increases from 2.5 x 10"'-* to 10"". 

The most striking change in behaviour occurs between the top 
panel and the second panel. At the shortest viscous timescale, the 
disc does not get into a trapped state, but evolves directly through 
corotation, while at lower viscosty the it always settles into the 
trapped state, with the inner edge moving in step with with the 
corotation radius. 

In all the discs, the initial evolution is the same: as m de- 
creases, Tin moves outwards, crossing rc over about 10^ years. 
Once that happens, however, the subsequent evolution of the disc 
and star differes substantially between simulations. In the most 
viscous discs (top), rin continues moving steadily outwards over 
10^" years - roughly 10^ times longer than the nominal spin-down 
timescale of the disc. Since rin keeps moving outwards, the torque 
on the star decreases too, so that the disc moves far away from rc 
before it is able to spin down the star 

As the viscosity is reduced, the ratio between the two 
timescales becomes smaller, and rin does not move so far away 
from rc before rin ~ rc. Thus for the disc with the lowest viscosity 
(bottom), rin and rc begin to move outwards after about 10^ years, 
and the trapped disc (where m is regulated by fc) has a much larger 
accretion rate onto the star (seven orders of magnitude larger after 
10* years) than for higher viscosities. 

By the initial condition chosen, the magnetic torque pushes 
the inner edge out across corotation. This causes mass to pile up 
outside rin. The higher the viscosity, the faster this pile is reduced 
again by outward spreading, and the faster the inner edge can con- 
tinue to move outward in response. The experimental result is then 



that trapping behavior is avoided when the transition takes place 
fast enough. We return to this in the discussion. 

The pile up is also influenced by the way in which accretion 
on the star changes as rin crosses rc, hence we expect that the pa- 
rameter controlling this, Ar2, will have a strong effect as well. 

In fig. |5] we show three discs in which the initial ratio 
TsD /Tvisc is kept fixed, but the value of Ar2 changes, from top to 
bottom, Ar2/rin = [0.4, 0.04, 0.004]. The larger Ar2, the further 
away rin must move from rc in order for the accretion rate to de- 
crease sufficiently to form the trapped disc. In the top panel (when 
Ar2 is largest), the disc must move out a considerable distance be- 
fore becoming trapped, sufficiently far to significantly decrease the 
efficiency of the spin-down torque (and hence increase the spin- 
down timescale of the star). As well as decreasing the spin-down 
efficiency, if Ar2 is large enough, the inner radius could move far 
enough away from rc that material could begin to be launched from 
the disc in an outflow, and the disc would become untrapped. 

As was shown in Section[3] the ratio TvIsc/Tsd itself can vary 
over many orders of magnitude in different systems, from 10~^^ in 
neutron stars with weak magnetic fields to 10~^ in discs around 
massive young stars. The size of this ratio will also determine 
whether a trapped disc can form. 

This analysis would suggest that, assuming Ar2/rin does not 
vary much from system to system, trapped discs are much more 
likely to form in protostellar discs than in strongly ionized discs 
around neutron stars. Furthermore, the closer rin is to rc, the higher 
the accretion rate in the trapped disc disc. Conceivably, especially 
if the viscosity in the disc were very low, this accretion rate could 
be larger than the average accretion rate in the disc itself, so that the 
disc could spin down the star for a long time without ever reaching 
spin equilibrium, even with a finite accretion rate onto the star 

As the results reported above show, the evolution can end ei- 
ther in a dead disc state in which the star has lost only a fraction 
of its angular momentum, or a trapped state in which corotation is 
maintained and the star can spin down much further Which out- 
come results depends details of the interaction between the star and 
the disc, parametrized in our model by the transition widths Ar and 
Ar2. It also depends on the rate at which the disc can respond vis- 
cously compared to the spin change rate of the star A fast response 
of the disc makes the transition through corotation faster than rc 
changes, and the disc is more likely to enter the dead state. This 
makes it far more likely to occur in young stellar systems than in 
neutron star binaries. In the results presented above, the initial con- 
ditions were taken from steady solutions of the viscous thin disc 
equation. These included dead discs in which a steady state was 
made possible by a sink of angular momentum at the outer bound- 
ary of the numerical grid, which takes up the angular momentum 
added by the magnetic torques at the inner edge. 

5.2 Dead discs evolving into trapped discs 

In this example we investigate the opposite case of section 15.11 
discs with rin initially outside corotation, and conditions chosen 
such that the disc begins by spreading inward. The mass flux at the 
outer boundary is set to zero. Varying the initial outer radius, rout,o 
varies the amount of mass in it, and the timescale of its long-term 
evolution can change. 

We adopt the representative model parameters used before, 
with an initial inner radius set to 1.3rc (corresponding to a neg- 
ligible accretion rate for our chosen Ar2/rin), and set rout,o ~ 
[10, 100, 10^, 10*] rin,o. The evolution of rin and rc is plotted in 
Figs.|6]and|7] Fig.[6]compares the evolution of rin/rc for different 
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Figure 4. The evolution from an accreting to a non-accreting disc, for increasing (top to bottom) ratios of T^i^c/TsT) (with Ar/rin = 0.1, Ar2/rin = 0.04). 
From top to bottom, the ratio Tvisc/TgD is 2.5 X [10"", 10~^, 10~^, 10~^]. As the ratio between the two timescales decreases, the disc is not able to move 
outwards as quickly before the star begins to spin down, so that r^^ will always remain close to r^. Black solid curve: evolution of r-^^. Red dashed curve: 
evolution of . 



sizes of disc, while fig.|7]compares the evolution of rin and rc for 
different disc sizes to the simulation where rout.o = 10rin,o. In 
both simulations, the different curves correspond to different initial 
Tout: 10 rin(dotted curve), 100 rin (dashed curve), 1000 rin (dash- 
dotted curve), lO^rin.o (dash-triple-dotted curve). 

Fig. |6] shows the evolution of r^-^jr^ in time for different ini- 
tial Tout- For the simulations with rout,o/'"in,o ~ [10, 100, 1000], 
the ratio ri^/rc declines to a minimum value that decreases as rout 
is taken larger, before again increasing approximately logarithmi- 
cally. In the largest disc, rout,o/'"in,o = 10'*, the evolution is the 
same as in the smaller discs at early times, but rin/r-c continues to 
decline for much longer until it reaches a minimum at around lO'tt 
when it finally turns over. 

The minimum value of rin/rc is thus determined by the 
amount of mass in the disc available for accretion. The larger discs 
have more mass, which sustains the accretion rate onto the star for 
a longer time before the drop in surface density causes rin to move 
outward again. 

The bottom panel of Fig. |7] shows the evolution of rin and 
Vc for rout.o = 10 rin,o. The evolution is qualitatively the same 
as we derived in the analytic approximation in Sec. 14.31 Initially 
Tin remains fixed as rc starts to evolve outwards, until moves 
close enough to rin that accretion can begin. The inner radius then 
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Figure 6. Discs starting outside corotation and spreading inward. Evolution 
of rin/rc for different initialrout. From top to bottom: rout,o = 10rin,o 
(dotted curve), 10^ (dashed curve), 10'' (dash-dotted curve), and 10* (triple 
dash-dotted curve) rj^ at t = 0. Larger discs have a larger reservoir of 
matter, so that they can sustain a larger m (so smaller rin) as rc increases 
due to spindown of the star. 
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Figure 5. The evolution from an accreting to a non-accreting disc, for stable discs (Ar/rin = 1) with different Ar2. From top to bottom, Ar2/rin = 
[0.4, 0.04, 0.004]. Curves show r^^ and Tc as in fig.|4] 



evolves outward at approximately the same rate as as the star 
spins down. The variation between different simulations is empha- 
sized in the top panel of Fig. |7] Here we plot the evolution of 
rc (thick curves, red) and rin (thin curves, black) for the discs 
with rout,o/nn,o ~ [10^, 10^, 10*], divided by the solution for 
''out,o/''iii,o = 10. For larger discs the accretion rate is higher, so 
that rin can move closer to r^. In the three smaller discs, the ac- 
creted mass adds a negligible amount of angular momentum to the 
star, so that as rin moves closer to rc. As a result, the spin-down 
torque simply becomes more efficient, and rc spins down faster. 
After IQ^t, for the disc with rout,o/rin,o = 10^ is more than 
10% larger than in the smallest disc. However, for the largest disc, 
the accretion rate puts rin close enough to rc that the spin-down 
torque starts to drop in efficiency and the spin-up from accretion 
becomes non-negligible. Although rc still increases, after lO'^t* rc 
is 30% smaller than for a small disc. 

These results show that size of the disc can considerably in- 
fluence the efficiency of spin-down, emphasizing the fact that the 
spin-down of a star is an initial value problem. The initial size of 
the disc can be as important as the ratio Tsd /Tvisc and the param- 
eter Ar2 in determining whether the disc can become trapped, and 
the efficiency of the spin-down torque. The results of this section 
would suggest that larger discs (with their larger reservoirs of mass) 
are more likely to become trapped than smaller discs. 



5.3 Long-term behaviour of discs of finite size 

As we found in section [5] the long-term evolution of the disc-l-star 
system tends to 'bifurcate'. The end state is either a star that has 
lost little angular momentum, surrounded by a dead disc, or a star 
continuously spun down by a disc trapped at corotation with the 
star. 

We investigate this further by a set of simulations in which 
the initial state is a disc of finite size ro, evolving in a grid that 
is 10 times larger, rout= 10 ro . Apart from this change we use 
the standard parameters (section 13.1) , setting the initial value of 
Tm/rc at 1.3 (for comparison with the results of the previous sec- 
tion). Whether the disc will be able to substantially spin down the 
star depends on the ratio the timescales Tvisc at ro to Tsd- If the 
disc is very small ( T^iBc(ro) <C Tsd), rin will move outward too 
quickly to spin down the star, while if rvisc(ro) S> Tsd, the mo- 
ment of inertia in the disc is sufficiently high to be able to absorb 
much more angular momentum from the star, so that the disc can 
operate as an efficient sink for the star's angular momentum. 

Fig. [8] shows results for three different sizes of disc. In each 
curve, the dashed red line shows the evolution of rc, while the 
solid black curve shows the evolution of rin. The top panel of fig. 
[8] shows the disc evolution when the viscosity is chosen such that 
7visc(7'o) -C Tsd, with ro = 100 nn. Since the angular momen- 
tum injected at rin is carried away by viscous spreading of the disc. 
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Figure 7. Evolution of rin and Tc in a dead disc witli different Tout ■ Bot- 
tom: Tlie evolution of rj^ (solid black curve) and r'c (dashed red curve) for 
the smallest disc, Tontfi = 10 fin.o- Top: The evolution of (thin black 
curves) and (thick red curves) for different sizes of disc, divided by the 
^out.o = 10 ?'in,o solution. The individual curves are the same as in Fig.|6] 



ally, however, viscous spreading of the disc wins, and rin begins to 
move outwards. 

The bottom panel of[8]shows the evolution of rin and Tc when 
T'visc(''o) 3> TsD- The result is essentially the same as in Sec. 
15.21 since the disc is now so large that additional angular momen- 
tum from the star does not reach rout on the spindown timescale. 
In other words, the moment of inertia of the disc itself is large 
enough that the star is able to spin down without causing rin to 
move rapidly outward. 

5.3.1 The rotation of Ap stars and magnetic white dwarfs 

An intriguing clue to the spindown of magnetic stars comes from 
slowly rotating Ap stars. As a class these stars are observed to 
have very strong dipolar magnetic fields (up to 10 kG). A few of 
them have extremely long rotation periods (up to 10-100 years), 
while others have rotation periods as short as 0.5 d. A similar phe- 
nomenon is observed in the magnetic white dwarfs. Most of these 
have periods of a few days or weeks, but some rotate as slowly as 
the slowest Ap stars. We suggest that the bifurcation of outcomes 
we found in the above is the underlying reason for the remarkable 
range of spin periods of magnetic stars. 
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Figure 8. Top: Evolution of Tin and r'c in a dead disc with mass transport 
through r'out, where Tvisc(''out,o) 2> Tsd- The moment of inertia in the 
disc is sufficiently large to prevent from moving out as Tc evolves. Mid- 
dle: Evolution of rin and in a dead disc with mass transport through 
i"out, where T'visc(''out,o) ^ TsT>- The disc is initially massive enough to 
spin down the star, but after some time evolves away from the equilibrium 
solution and moves rapidly out. Bottom: Evolution of rjn and in a 
dead disc with mass transport through rout, where Tvisc (rout,o) ^ Tsd- 
The moment of inertia in the disc is sufficiently large to prevent rj^ from 
moving out as rc evolves. Curves show rjn and rc as in fig.O 

the disc quickly evolves away from its initial configuration. Since 
the moment of inertia in the disc is much smaller than in the star, 
the disc becomes too spread out to absorb the angular momentum 
injected at rin, and rin moves outward before the star is able to slow 
down substantially. 

In the middle panel of [S] we show the evolution of rin and 
rc when the spin-down timescale is comparable to the viscous 
timescale at ro. The disc also initially diffuses outwards (as seen 
by the increase in rin around i = 5 x 10"'), so that the rate of angu- 
lar momentum exchange from the star to rin decreases. This allows 
rc to catch up, and the two radii start to evolve together Eventu- 



6 CONCLUSIONS 

As found before in lSunvaev & Shakural ( Il977h and DSIO, a disc in 
contact with the magnetosphere of a star can be in a 'dead' state, 
with its inner edge well outside the corotation radius so the ac- 
cretion rate onto the star vanishes, and the torque exerted by the 
magnetic field transmitted outward by viscous stress. In the calcu- 
lations reported in DSIO, an additional state was found, interme- 
diate between the accreting and dead state. In this state, the inner 
edge of the disc stayed close to the corotation radius rc, even as 
the accretion rate onto the star varied by large factors. We call this 
phenomenon trapping of the disc. Accretion can be stationary or in 
the form of a limit cycle in this state. 

One of the goals of this investigation was to find out under 
what conditions this trapping takes place, using a series of numer- 
ical experiments with varying initial and boundary conditions. If 
initial conditions are such that the inner edge starts inside corota- 
tion and slowly moves outward, we find that the disc gets stuck 
in a trapped state for a long time if the disc viscosity v is low 
(up to about 10'^ times shorter than the spindown timescale of the 
star). The accretion rate then slowly vanishes but the inner edge 
always stays close to corotation. At higher viscosity (> 10^ x the 
spindown timescale) on the other hand, the disc evolves through 
corotation into a 'dead' disc state, with inner edge well outside the 
corotation radius. In terms of the standard viscosity parametriza- 
tion V = a{H/r)^, compared with the spin-down timescale of the 
disc, a trapped state is more likely to happen in strongly magnetic 
(Bs ~ 10^^ G) X-ray pulsars and protostellar discs than the weaker 
magnetic fields (Bs ~ 10** G) of millisecond X-ray pulsars. 

Our second goal was to find out how a star spins down in the 
long term, under the influence of the angular momentum it loses 
by the magnetic torque exerted on the disc. The results show an 
interesting 'bifurcation' of long-term outcomes: if the disc evolves 
into a dead state, the star loses only a fraction of its initial angular 
momentum, and can remain spinning rapidly throughout its life. If 
on the other hand it enters a trapped state at some point, it remains 
in this state. The star can then slow down to very low rotation rates, 
the inner edge of the disc tracking the corotation radius outward. 
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We suggest that these two outcomes can be identified respectively 
with the rapidly rotating and slowly rotating classes of magnetic 
Ap stars and magnetic white dwarfs. The evolution of the trapped 
state could also be reproduced with a simplified model that does 
not require solving the full viscous diffusion equation. 

This picture of magnetosphere-disc interaction differs from 
the standard view that mass will be 'propellered' out of the sys- 
tem instead of accreting, once the star rotates more rapidly than the 
inner edge of the disc. As shown already b y Sunvaev & Shakural 
( Il977h and again argued in lSpruit & TaamI Il993h and DSIO this 
assumption is not necessary, in many cases unlikely, and ignores 
some of the theoretically and observationally most interesting as- 
pects of the disc-magnetosphere interaction. 

Though the dead stateis a regime where a significant fraction 
of the disc mass could in principle be expelled (rin2> ^c), the re- 
sults presented here show that magnetospherically accreting sys- 
tems often avoid this regime. Instead, they end up in the trapped 
state, in which the disc-field interaction keeps the inner radius trun- 
cated very close to the co-rotation radius, even at very low accretion 
rates. Both the trapped and (rin~ r^) and the dead state (rinS> r^) 
allow the disc to efficiently spin down the star. The disc retains a 
large amount of mass, but in the absence of accretion onto the cen- 
tral star appears quiescent. 



7 DISCUSSION 

By assuming a given dependence of viscosity on distance, we have 
bypassed the physics that determines it. In terms of the standard a- 
parametrization of viscosity, we have left out the physics that deter- 
mines the disc temperature and hence its thickness H/r. Additional 
time dependence or instabilities may arise from feedback between 
accretion and disc temperature. The radiation produced by matter 
accreting on the star could be large enough to change temperature, 
the ionization state and hence the thickness at the inner edge of the 
disc. In a trapped disc, this is just the region that controls the accre- 
tion rate onto the star. The size of the transition region, which we 
have parametrized with the widths Ar, Ar2 may well depend on 
disc thickness. Positive feedback may be thus possible. We leave 
this possibility for future work. 

In systems such as the accreting X-ray pulsars the accretion 
is episodic on long timescales. This is attributed to the instability 
of viscous discs that is also responsible for the outbursts of cata- 
clysmic variables. In these cases, in which there is a large change 
in rh and rin is far from in the quiescent phase, there is the pos- 
sibility of hysteresis: the same accretion rate will lead to a different 
value of Tin (and therefore disc torque) depending on whether the 
source is moving into or out of outburst. As the source goes into 
outburst, the disc will not have as much mass in its inner regions, so 
that Tin will move inward gradually from large radii until it crosses 
Tc and the outburst begins. In the decline phase, however, the disc 
will become trapped around the inner radius of the disc when the 
accretion rate drops, allowing for a larger spin-down in the disc and 
accretion bursts via the instability of DSIO. The net effect of such 
episodic accretion on the spin history of the star, as compared with 
the case of steady accretion, is not obvious. We discuss this in more 
depth in the companion paper. 

Some work on disc-magnetosphere interaction assumes mag- 
netic torques to act over a significant part of the disc(K6nigl 1991 ; 
lArmitage & Clarke! Il996l) . More recent work (and in particular 
numerical MHD simulations of the disc-field interaction) finds 
the interaction region to be much narrower (as we have also as- 



sumed here), and spindown torques on the star are correspondingly 
smaller. The difference for the long-term spin evolution is not dra- 
matic, however, as our trapped disc results demonstrate. 

Interestingly, lArmitage & Clarkel d 1 996b also observed that 
their discs would become trapped around as the accretion rate in 
the disc decreased by several orders of magnitude. In their model, 
the magnetic field-disc interaction also acts like a boundary con- 
dition at low accretion rates and the disc evolved viscously in re- 
sponse. This is presumably because although the disc in their model 
is threaded by a magnetic field everywhere, the disc-field interac- 
tion is by far the strongest in the inner regions, so that they see a 
similar behaviour to the one described in this paper. 

MHD simulations of interaction between the disc and the mag- 
netic field are becoming increasingly realistic. These simulations 
can only run for very short timescales (the longest of order Tvisc at 
Tin), so they tend to emphasize initial transients. Still, they offer in- 
sight in how the disc and magnetic field will interact. To date, most 
simulations have concentrated on strongly accretin g dHavashi et al.l 
'1996; Goodson et al. 1997'; iMiller & Sto"n3ll997h or propellering 
(Romanova et al. 2004) cases. Simulations that come closest to the 
conditions of a tr apped disc are the study of a so-called 'weak pro- 
peller' regime by lUstvugova et aP (|2006|) . The authors found that 
discs in which rin was initially truncated close to launch much 
weaker outflows than discs truncated further away. They also found 
some evidence of the field changing the disc structure (as shown 
from mass piling up in the inner regions), although in their simula- 
tion the majority of the angular momentum in the star was carried 
away via a wind, rather than through the disc. These simulations 
also did not run for very long, however, so it is hard to separate 
transient behaviour due to the initial conditions from the longer- 
term systematic effects we are interested in. 

iPema et al. did calculations were done for a case where 

the star's field is inclined with respect to the disc axis. The au- 
thors' results suggest that the transition width which we have 
parametrized by Ar2 would increase with the inclination of the 
magnetic field, since there will be some values for rin at which 
which both accretion and disc mass trapping can occur This then 
would suggest that dead discs are more likely to form in sys- 
tems with small inclinations between the spin axis and magnetic 
axis, and also that the disc instability studied in DSIO (which also 
tends to occur for smaller values of Ar2). This prediction is sup- 
ported by the observation that Ap stars with the longest periods 
tend to have the lowest inc lination angles for the magnetic field 
jLandstreet & MathvslbOOOl) . 

The transition widths of disc-magnetosphere interaction that 
can be inferred from these simulations are significant, and are in the 
range we have assumed here. They are much larg er than the very 
narrow interaction regions assumed by Matt et al. dMatt & Pudritj 
2004, 2005; Matt et al. 2010). 

We have found that the distance from at which the inner 
edge of the disc gets trapped is determined by the ratio of two 
timescales for the disc's evolution: the viscous evolution timescale 
in the disc (which determines the rate at which disc density profile 
can change) and the star's spin-down timescale (which sets the rate 
at which moves outwards). The viscous timescale is in general 
much shorter than the spin-down timescale, and the ratio of the two 
timescales varies from ~ 10^'' (in protostellar discs) to ~ 10^^^ 
(in millisecond X-ray pulsars). The larger the ratio, the longer the 
disc will take to respond to changes in Tc. As a result rin remains 
closer to Vc than it would if the ratio were smaller, and the trapped 
disc has a higher accretion rate onto the star. If the ratio is too low. 
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then Tin moves much further away from r^, and the system likely 
enters the dead disc regime. 

For the parameters characterizing the interaction region be- 
tween disc and magnetosphere in our model, we have used here 
values such that the cyclic accretion behavior found in DSIO does 
not develop. This was done for convenience, since the short time 
steps needed to follow these cycles makes it harder to calculate the 
long-term evolution. The time-averaged effect of these cycles is not 
expected to make a big difference for the long-term evolution. 

These cyclic accretion bursts can persist in the trapped state, 
when the star is spinning down efficiently. They are observed to 
occur both over several orders of magnitude of accretion rates and 
transiently (over a small range of accretion rate), and depend on 
the ratio of timescales discussed above. Whether or not the insta- 
bility occurs is determined by the detailed disc-field interaction 
(the parameters Ar and Ar2 in our model). The peak of the ac- 
cretion bursts is typically much larger (> 10 x) than the mean 
accretion rate for the system, and the period is typically between 
0.01 — 10'^rvisc(''in). The properties and conditions for occurrence 
of these cycles are studied further in a companion paper 



7.1 'Propellering' 

In our calculations we have left out the possibility that interaction 
of the magnetosphere with the disc will cause of mass ejection from 
the system. The point being that, contrary to common belief, such 
interaction can functio n without mass ejection by ' propellering', as 
pointed out already bv lSunvaev & Shakural digTVh . Understanding 
of this restricted case, as we have developed here, is prerequisite 
for understanding the much less well defined case of mass loosing 
discs. 

On energetic grounds, mass loss from the system is nec- 
essa rily limited, unless th e inner edge is well outside corota- 
tion jSpruit & Taanil 1 1993 1). This point has also been made by 
iPema et al. I ( l2006h . who propose that mass lifted at rin may fall 
back on the disc at some finite distance. This would create a feed- 
back loop in the mass flux through the disc, opening the possibility 
of additional foms of time-dependent behaviour that do not exist 
in accretion onto non-magnetic stars. In the trapped disc state we 
have studied here the difference in rotation between disc and star 
is small, so any significant amount of mass kicked up from the in- 
teraction region cannot move very far before returning to the disc. 
Its effects are then secondary, at least for the long-term evolution 
of the disc. 

The possibility of significant effects of mass loss is more re- 
alistic for the dead disc states, where the distance of the inner edge 
from corotation can become much larger. 

Real propellering is expected to happen when mass transfer 
from a companion star sets in for the first time onto a rapidly spin- 
ning magne tic star. The catacly smic binary AE Aqr is evidently 
such a case jPearson et al.ll2003l) . A disc is absent in this CV, and 
all mass transfered appears to be ejected in a complex outflow as- 
sociated with strong radio emission. 
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